Kohn anomalies and non-adiabaticity in doped carbon nanotubes 
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The high-frequency Raman-active phonon modes of metallic single-walled carbon nanotubes (SWNTs) are 
thought to be characterized by Kohn anomalies (KAs) resulting from the combination of SWNTs intrinsic one- 
dimensional nature and a significant electron-phonon coupling (EPC). KAs are expected to be modified by the 
doping-induced tuning of the Fermi energy level ep, obtained through the intercalation of SWNTs with alkali 
atoms or by the application of a gate potential. We present a Density-Functional Theory (DFT) study of the 
phonon properties of a (9,9) metallic SWNT as a function of electronic doping. For such study, we use, as 
in standard DFT calculations of vibrational properties, the Born-Oppenheimer (BO) approximation. We also 
develop an analytical model capable of reproducing and interpreting our DFT results. Both DFT calculations 
and this model predict, for increasing doping levels, a series of EPC-induced KAs in the vibrational mode 
parallel to the tube axis at the T point of the Brillouin zone, usually indicated in Raman spectroscopy as the G~ 
peak. Such KAs would arise each time a new conduction band is populated. However, we show that they are an 
artifact of the BO approximation. The inclusion of non-adiabatic (NA) effects dramatically affects the results, 
predicting KAs at T only when 6f is close to a band crossing Ex ■ For each band crossing a double KA occurs 
for 6f ~ Ex ± hto /2, where hui is the phonon energy. In particular, for a 1 .2 nm metallic nanotube, we predict 
a KA to occur in the so-called G~ peak at a doping level of about N e i/C = ±0.0015 atom (ef ~ ±0.1 eV) 
and, possibly, close to the saturation doping level (iV e ;/C~ 0.125), where an interlayer band crosses the n* 
nanotube bands. Furthermore, we predict that the Raman linewidth of the G~ peak significantly decreases for 
\cf\ > hu/2. Thus our results provide a tool to determine experimentally the doping level from the value of 
the KA-induced frequency shift and from the linewidth of the G~ peak. Finally, we predict KAs to occur in 
phonons with finite momentum q not only in proximity of a band crossing, but also each time a new band is 
populated. Such KAs should be observable in the double-resonant Raman peaks, such as the defect-activated D 
and peak, and the second-order peaks 2D and 2G. 



> 

m 

On 

(N 
O 

o 

-i— > 



i 

C 

O 

o 



X 
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I. INTRODUCTION 

Since their discovery in 1991—, carbon nanotubes have 
raised an enormous interest both from the academic and the 
technological points of view. They exhibit in fact a variety of 
exciting features: their quasi-one-dimensional nature, due to a 
diameter of 1-2 nm and a length of up to several micrometers, 
makes them sharp probes for scanning tunneling microscopes 
and an excellent model for one-dimensional physics. Their 
mechanical and tensile strength make them of great interest 
in composite materials, and, most significantly of all, they 
have very unusual and extremely promising electronic proper- 
ties, displaying metallic or semiconducting behavior accord- 
ing to their structure and helicity^. The electronic proper- 
ties of SWNTs, already particularly interesting from the tech- 
nological point of view, promise to be the future of nano- 
electronics due to their tunability, achieved by doping the nan- 
otubes through intercalation with alkali atom s 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 ' 12 ' 13 
or application of a gate potentia l 14 ' 15 ' 16 ' 17 . However, there are 
still a number of experimental challenges to be solved in order 
to fully develop SWNT-based nano-electronic technology, in 
particular the low-cost industrial-scale synthesis of nanotubes 
of good chemical purity, crystalline quality and given helicity. 

An experimental tool largely used to characterize SWNTs is 
Raman spectroscop y 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 15 ' 16 ' 17 ' 18 ' 19 ' 20 ' 21 ' 22 ' 23 . Typical Ra- 
man spectra of carbon nanotubes display a peak around 150- 
300 cm -1 , due to the radial breathing mode (RBM), which 



has been recently used as a tool to infer the size and chiral- 
ity of the nanotubes 20,23 . Other important features of SWNT 
Raman spectra include a peak around 1350 cmT 1 , activated 
by defects and impurities, and known in the literature as D 
peak, and a large structure around 1570 cm -1 , due to modes 
tangential to the nanotube and known as G peak. This last 
feature is thought to have two components, usually referred 
to as G + and G~, originating from the E2 9 in-plane modes 
of graphite. In refs i 24 ' 25 ' 26 it has been shown that in metallic 
SWNTs the G + component corresponds to the tangential vi- 
brational mode perpendicular to the nanotube axial direction, 
while the G~ component corresponds to the tangential vibra- 
tional mode parallel to the nanotube axial direction. In the 
following, we will refer to the former as "the nanotube TO 
tangential mode", and to the latter as "the nanotube LO axial 
mode". Some Raman studie s 5 ' 6 ' 7 ' 8 show that the frequency of 
the G peak increases up to 1600 cmT 1 at low doping levels, 
and then suddenly drops to about 1550 cmT 1 at the saturation 
threshold of alkali intercalation, estimated around a number 
of electrons per carbon atom N e i/C — 0.12 (MCs). 

Other experimental 916,18 and theoretica l 27 ' 28 ' 29 works on the 
effect of doping on SWNTs report a similar G-peak soften- 
ing or even a Luttinger-Fermi liquid transition^ 2 -. Since the 
LO axial and TO tangential modes are particularly sensitive 
to the electronic structure of SWNTs around the Fermi energy 
24 ' 26 i 27 ' 30 ' 31 ep, and ep directly depends on the charge doping 
level, a profound understanding of the interplay between the 



2 



vibrational and the electronic properties of nanotubes looks to 
be crucial for technological development. 

In this work we report our theoretical study of the elec- 
tronic and vibrational properties of doped SWNTs, based on 
DFT first-principles calculations and analytical results. We 
will show in the following that: i) the vibrational properties of 
SWNTs can be obtained from the so-called electronic zone- 
folding of a graphene sheet; ii) their behavior as a function 
of the charge doping level can be determined by the knowl- 
edge of electron-phonon coupling (EPC) in graphene; Hi) or- 
dinary quantum-mechanics calculations relying on the adia- 
batic Born-Oppenheimer approximation fail when applied to 
SWNTs, where non-adiabatic effects are enhanced by their 
intrinsic one-dimensional nature. 

Our paper is organized as follows: in section HI] we will 
describe our theoretical framework, and in particular how the 
Raman-active modes of a nanotube can be accurately obtained 
through an appropriate electronic sampling of the graphene 
Brillouin Zone (BZ). We will then show in section [III] that 
the DFT results can be almost perfectly reproduced by an in- 
tegral model, that uses the graphite EPC and the slope of the 
electronic bands around the (undoped) Fermi level e F (usually 
referred to as the ir and it* bands) as external inputs, and that 
becomes analytical in the limit of vanishing temperature. The 
results on the DFT and the model-derived LO axial and TO 
tangential modes of a SWNT will be presented, showing that 
Kohn anomalies (KA) would occur, within the adiabatic ap- 
proximation, each time a new electronic band is populated by 
the electrons at ef. In section|IV]we show that when the Born- 
Oppenheimer approximation is lifted the outcome is dramat- 
ically different, and that a drop in the frequency of the G~ 
peak occurs at a doping level such that e F — e F « 0.1 eV or 
close to the saturation level. Section [V] will be devoted to the 
discussion of the physical properties of metallic SWNTs that 
are experimentally accessible, and to the conclusions. 

II. THEORETICAL BACKGROUND 
A. General description of graphene and SWNTs 

Graphene is a semimetal, and its highest valence and lowest 
conduction bands have a circular conical surface shape around 
the points K and K' of the hexagonal first Brillouin zone at 
zero doping; its Fermi surface is reduced to the vertexes of the 
cones at those two points. 

The structure of a SWNT is obtained by rolling a sheet of 
2D graphene, and closing both ends with fullerene-like semi- 
spheres. However, nanotubes are commonly studied as in- 
finite ID crystals, i.e. neglecting the details of those ends. 
Nanotubes are then uniquely determined by the knowledge of 
the chiral vector Ch = na^ + ma.?,, where ai and a? are the 
basis vectors of graphene, n and m are integers 2,3 . The trans- 
lational vector determines the periodicity along the tube axis 
and is T = ( 2 "?+") ai _ ( 2 "+ m ) a2, where d^ is the great- 
est common divisor of (2m + n) and (2n + to). Nanotubes 
with indexes such that n — to = 3fc, where k is an integer, are 
metallic, otherwise they are semiconductors. Nanotubes with 



identical indexes, i.e. n=m, are thus always metallic, and are 
usually called "armchair" tube. 

Electronic zone folding (EZF) is a common approximation 
for the electronic band structure of a SWNT, which consists 
in neglecting the curvature of the SWNT. Within the EZF the 
electronic states of a SWNT are approximated by the elec- 
tronic states of a graphene sheet having wavevector k such 
that k • Ch = 2-kv, being v an integer. The wave-vectors 
allowed for the SWNT satisfy 

k = i/k ± + Ak||/||k|||| (1) 

where k is a real, k|| is parallel to the tube axis, k^ is per- 
pendicular to the tube axis and satisfies ||kj_|| = 27r/||Ch|| — 
2/d, d being the tube diameter. 

This way of constructing the electronic band structure ac- 
counts for the fact that SWNTs can be either metallic or semi- 
conducting. In fact, if in SWNTs a line of allowed wave vec- 
tors k crosses the K or the K' point, the nanotube is metal- 
lic, otherwise it is a semiconductor. In particular, in the BZ 
of armchair metallic SWNTs the section of the vertical plane 
containing the v = line of k vectors and the conical bands 
consists of two straight lines of slope ±/3, while the conic sec- 
tions for \u\ > 1 consist of arms of hyperbolae. The energy 
difference between the maximum of the lower arm and the 
minimum of the upper arm are usually indicated as E^ v in the 
metallic case (see also Fig.O. 

Previous work s 2,3 i 24 - 25 ' 32 have shown that the effects of the 
nanotube curvature can be neglected as a first approximation 
for SWNTs whose diameter is larger than 1 nm (the most 
common nanotubes present in experimental samples) and the 
electronic properties of a nanotube can be quite accurately ob- 
tained from EZF. 



B. Electronic zone-folding for phonon calculations 

The explicit calculation of the full dynamical matrix of an 
infinite nanotube, even within an efficient scheme such as 
Density-Functional Perturbation Theory (DFPT) 33 , is a quite 
demanding task from the computational point of view. On the 
other hand, the SWNT vibrational modes typically observed 
in Raman spectra can be traced back to the phonon and elastic 
modes of an isolated graphene sheet. 

However, it has been shown that the phonon dispersions 
of SWNTs cannot be as accurately reconstructed from the 
graphene ones as the electronic bands, because of the presence 
of KAs which behave differently in graphene 34 and metallic 
SWNTs 24,25,26.27,30.31^ The KAs &xe determined by singular- 
ities in the electron screening and such singularities depends 
on the dimensionality of the electron Brillouin zone. The di- 
mensionality is two for graphene and one for SWNTs. This re- 
duced SWNT dimensionality is due to the quantization (con- 
finement) of the electronic wavevector around the nanotube 
circumference, which is actually described by the EZF. This 
last consideration suggests a practical scheme for the calcula- 
tion of the phonon dispersion in SWNTs which neglects cur- 
vature, but fully takes into account the more important con- 
finement effects. Phonons in nanotubes are well approximated 
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by the phonons of a fiat graphene sheet, if the calculation is 
done performing the electronic Brillouin-zone integration on 
the lines of the electronic zone-folding (Eq.[TJ. Such phonon 
calculation method (phonon-EZF) was introduced in refs i 24 ' 25 . 
A phonon-EZF calculation requires the use of a unit cell con- 
taining two atoms and thus, is clearly much less computa- 
tionally demanding than full calculations on an actual SWNT, 
which requires a unit-cell with tens or hundreds of atoms. In 
Subsect. HID I we will demonstrate that, for the preset study, 
curvature effects are negligible and the phonon-EZF provides 
an accurate description of the high-frequency phonon modes. 



C. Computational details 

Our first-principles calculations are based on DFT, within 
the plane-wave (PW)/pseudopotential scheme implemented in 
the Quantum-ESPRESSO code 35 . We adopt a Perdew-Burke- 
Erzherhof gradient corrected functional, and an ultrasoft pseu- 
dopotential to describe the C atom. A PW kinetic energy cut- 
off of 30 Ry is sufficient to ensure convergence on the struc- 
tural, electronic, and vibrational properties. We choose as a 
case-study the (9,9) metallic armchair SWNT, containing 36 
distinct carbon atoms, and having a tube diameter of 1 .24 nm. 
The full nanotube calculations are performed by setting the 
nanotube in an infinite lattice of hexagonal symmetry in the 
plane perpendicular to the tube axis, mimicking thus the bulk 
effect of a real nanotube bundle or an isolated tube by tun- 
ing the lattice parameter between 1 .56 nm and 1 .72 nm. The 
c axis is maintained constant. Integrations in the nanotube 
one-dimensional Brillouin zone have been performed by using 
regular grids of k points along the c reciprocal axis. We use a 
Fermi-Dirac electronic smearing of 0.01 Ry (corresponding to 
an electronic temperature of 1578 K) for SWNT full phonon 
calculations, and of 0.002 Ry (315 K) for structural calcula- 
tions and frozen-phonon tests (see below). Grids of 8 and 40 
k points were sufficient, respectively, to ensure a good con- 
vergence (~ 5 cm -1 ) on the LO axial frequencies in the two 
cases. 

In the phonon-EZF calculation we use a graphene hexago- 
nal unit cell of lattice parameter a cxp = 2.47 A and c = 5 A. 
The number of points is chosen as to ensure phonon frequen- 
cies converged within about 0.3 cm -1 , and is inversely pro- 
portional to temperature. At a temperature of 315 K, the total 
number of points is 3240, equivalent to 180 points along the c 
direction of the reciprocal space of the (9,9) nanotube, while at 
1578 K this number reduces to 720, equivalent to 40 points in 
the nanotube case. For both the full nanotube and the phonon- 
EZF case, the effect of (low) charge doping is simulated by 
adding an excess electronic charge which is compensated by 
a uniformly charged backgournd. This is done using the stan- 
dard implementation of the ESPRESSO-package 35 . This ap- 
proach is known to describe very well the properties of doped 
SWNTs 36 . This is justified even in the case of alkali interca- 
lation by the report that at low doping levels, i.e. up to MC15, 
the electronic charge of the metal atoms is thought to be com- 
pletely transferred 37 , and that the effect of the intercalation is 
essentially described by taking into account the sole effect of 
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FIG. 1: Variation of the LO axial frequencies with electron doping at 
315 K, calculated from phonon-EZF (full diamonds) or through full- 
nanotube frozen-phonon (open circles). The solid line is a guide to 
the eye. The metal/carbon composition corresponding to the doping 
level is indicated on the top of the figure. Inset: variation with elec- 
tronic doping of the phonon-EZF LO axial frequency (full diamonds) 
compared to a SWNT full DFPT dynamical matrix calculation (open 
squares), at 1578 K. 



the charge transfer. 



D. Accuracy of the method of phonon calculation 

The first step of our study is to assess the accuracy of the 
phonon-EZF method with respect to the SWNT phonon cal- 
culations. In Fig. [T] we show the frequency of the LO axial 
mode as a function of electronic doping at 315 K, obtained 
both through the phonon-EZF method or by an explicit nan- 
otube frozen-phonon calculation. In the latter case, the atoms 
are displaced along the axial direction of the nanotube, and 
the frequency is obtained by the curvature of the displacement 
energy. The convergence of this single-mode frozen-phonon 
frequency has been confirmed by testing the result with a more 
accurate sampling of the BZ, that is up to 200 uniformly- 
spaced k points in the one-dimensional BZ grid. In the inset 
we analogously compare the phonon-EZF LO axial frequency 
with a full DFPT nanotube phonon calculation at 1578 K. In 
this latter case, all the 108 T vibrational modes of the nanotube 
are explicitly calculated; a lower temperature, and thus a cor- 
respondingly larger grid of at least 40 k points, would be com- 
putationally too demanding for a complete and accurate study 
of the problem. Both graphs show that phonon-EZF not only 
determines the nanotube phonon frequencies in good agree- 
ment with explicit SWNT calculations, but it also captures 
very well the strikingly non-monotonic dependence of the LO 
axial mode with respect to electronic doping, to be discussed 
later on, confirming that the effects of the SWNT curvature 
can be confidently neglected for a (9,9) armchair nanotube, 
and that phonon-EZF is a reliable and accurate method. 
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FIG. 2: Variation of the electronic doping per C atom as a function 
of the Fermi energy ep. Solid line: analytical expression (Eq|2j with 
En/2 = 0.83 eV; diamonds: DFT calculations, within phonon- 
EZF, of the Fermi energy in a (9,9) SWNT at 315 K. The E?[/2 and 
E22/2 levels are indicated by vertical dotted lines. Inset: Calculated 
electronic bands of the (9,9) SWNT around the K point, is taken 
as the zero energy, Ef( /2 and E22 /2 are the minima of the second 
and third (hyperbolic) conduction bands. The valence bands (below 
e%) and the conduction bands (above e%) are usually referred to as 
the 7T and n* bands, respectively. The bands crossing at e a F have the 
index v — 0, while the valence and conduction bands immediately 
below and above have indexes v = ±1 (see Eq.QJ. 



We report in Fig. [2] the corresponding curve relating ep 
to N e i/C, along with our DFT calculations at 315 K. The 
most significant discrepancy being between the ideal value of 
Eff/2 and the calculated one, which is 0.83 eV. If however 
we treat Ef[/2 as an independent parameter, and we red-shift 
it by 0.06 eV, the agreement is very good. This indicates 
that the analytical expressions can be safely used to develop 
an integral and analytical model, based on the EPC and the 
electronic bands, aimed at describing and predicting the be- 
havior of the LO axial and TO tangential modes in metallic 
nanotubes. The electronic bands of the nanotube, explicitly 
calculated by DFT, are shown in the inset around the K point 
of the BZ, along with the band nomenclature tt/tt* , v used in 
the literature and adopted in this work. One can notice that 
the minimum of the v = ±1 conduction band is slightly dis- 
placed to the right, and that the bands deviate from the conical 
shape far from K. In the following we will use a polynomial 
fit of the calculated bands in numerical evaluation of integral 
expressions, while the ideal conical shape of the bands will be 
used to determine the analytical limits at low temperature. 



III. INTEGRAL AND ANALYTICAL MODEL 

A. Electron-phonon coupling contribution to the dynamical 
matrix 

The frequency of a vibrational mode at a q wavevector, is 
obtained from the dynamical matrix D q 



E. Fermi level in doped SWNTs 

The key quantity intervening in the anomalous vibrational 
properties of SWNTs is the Fermi energy level ep. In un- 
doped graphene and in metallic armchair nanotubes, its value 
(e F ) coincides with the band energy at the crossing of the it 
and 7r* bands at the K and K' points of the graphene BZ, and 
is chosen as the zero band energy hereafter. Since the results 
are independent on the sign of ep, we consider here and in the 
following that ep > 0, without loss of generality. Following 
the description of the graphene electronic bands in terms of 
conic sections as in subsection III Al the increase of ep due to 
electronic doping per carbon atom N e i/C in a (n, n) SWNT 
can be analytically expressed, at low doping levels and a tem- 
perature T=0 K, by the following formula: 



N el /C(e F ) 



qfvg 

n(3d 



e F +9(e F -E?{/2) 



2aQ\/3 
' n(3d 



(3) 



where M is the atomic mass of carbon. As shown in refs. 25 ' 34 , 
the dynamical matrix can be written as the sum of an EPC 
direct contribution Z) q , which contains non-analytical terms 
giving rise to Kohn anomalies, and a term containing all the 
other contributions: 



D n = D r 



D, 



other 



(4) 



The direct contribution of EPC to the dynamical matrix can 
be written as 



q Mc ^ 



/(<*,») - /(fk+q,/) 
£k,i ~ e k+q,/ 



• |G(k+q),/;k,i| 2 (->) 



(2) 

where clq is the lattice parameter, j3 is the slope of the conical 
bands and equals 14.1 eV^, and 9(x) is the step function. As 
previously mentioned, the energies of electronic transitions in 
SWNTs are usually indicated as Eff^ S , where, for armchair 
SWNTs, i coicides with the v band index of the initial and 
final states, and the superscript refers to metallic or semicon- 
ducting SWNTs. In the case of a (n, n) nanotube, assuming 
perfect conical bands, E\{j2 = 2/3/ 'd. In particular, for the 
(9,9) nanotube Ef( /2 = 0.89 eV. 



-where is the number of points in the SWNT Brillouin 
zone, i and / are the band indexes indicating the two states 
involved in the the electronic transition. Since we consider 
only the contribution of the tt and ir* bands, i, f — ir, ir*. The 
function /(e) = {exp[(e— ep)/A:T] + l}~ 1 , is the Fermi-Dirac 
distribution; G(k +q ) ,/ : k,i is the EPC matrix element, defined 
as 



G( k+q) , /;k , t - (k + q,/|Ay q |k, 



(6) 



where AVq is the first-order derivative of the Kohn-Sham 
self-consistent potential with respect to the atomic displace- 
ments corresponding to a q-vector phonon; |k, i) is the Bloch 
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electronic wavefunction. Hereafter, we develop our model 
on the basis of the EZF sampling, and the k vectors are 
those allowed by Eq. Q] Since we study phonons close to T 
we limit ourselves to wave vectors entirely along kii writing 
q = <7k||/||k|| || and ||q|| <C ||K||. Therefore v is conserved 
through the transition i — > /. In the limit of a nanotube of in- 
finite length we replace ^ k by 2~^„ J dk, which yields 
for the non-analytical part of the dynamical matrix 



where we distinguish between interband (i ^ f, first line) and 
intraband (i = /, second line) transitions, for which we used 
the following limit: 



lim 



f[^f{k')]-f[e vJ {k' + q)] _ df 



Then if the temperature T — > 0, 



,,/(&')]• (13) 



T 



E 



dk 



f[e„4k)}-f{e„j(k + q)} 
e v ,i{k) - e v j{k + q) 



\G 



(k+q),/;k,i| 
(7) 

Neglecting the effects of curvature, the EPC matrix element 
(7(k+q),/ ; k,i of a nanotube of diameter d and longitudinal pe- 
riod T is related to the EPC G(k+q),/ ; k,j of graphene by the 
ratio of the unit cells areas (Eq.4 o^&): 



de 



e»Ak')] 



4n 2 d 
V3a 2 £ 



E 

± 



DOS v , 



d7rT|G( k+q ) j. ki . 



|G(k+q),/:k,i| 2 (8) 



The main contribution to the EPC term originates from the K 
point of the graphene BZ, so we can define the wavevector 
k' = k K. Furthermore, since KAs originate when the 
denominator in Eq. [7] vanishes, we can restrict the integral in 
the BZ to a small interval of width 2k around K (and K'). We 
can also define the angles 9, between k' and q, and 9', between 
k' + q and q. In refs. 24 ' 34 it was shown that for phonons close 
to r, the nanotube EPC can be in principle expressed in terms 
of the graphene EPC (G r ) F = 45.60 (eWA) 234 , modulated 
by a geometric factor: 



a {e F )5(k' -k F vJ J 
(14) 

where D0S^j i(7 (e F ) and k F f a are the electronic density of 
states (DOS) per G atom and the Fermi wavevector of the 
/ band of v index and slope of sign a. The intraband EPC 
contribution to the dynamical matrix at T is thus proportional 
to the density of states at the Fermi level. 

The dynamical matrix at finite q can be calculated by using 
the general expression (fTTT) : in the following we will focus on 
the contribution at q=0 of the lowest conduction band v=Q, 



£)p, and the second lowest bands v 
of the Fermi energy 



±1, -Dp 1 , as functions 



Dr(e F ) = DUe F ) + Dp 1 (e F 



(15) 



,^TO/LO ,2 _ i n 
^(k+q),/;^! - \ U 



) F [1 ± sign(e„j ■ e Vji ) cos(6» + 6')] 

(9) 



The contribution of higher conduction bands could be in prin- 
ciple included with no difficulty in the same form as the 
v = ±1 bands, but their effect is negligible for the doping 
levels considered in this work, up to N e i/C ~ 0.06, i.e. the 
value of estimated total charge transfer in the case of alkali 
intercalation 3 -^ 



where 



cos 



9') = 



k'(k' + q) - (yk^f 



y/k' 2 + (vk ± ) 2 ^(k> + q) 2 + (isk ± ) 2 



(10) 



and the +(— ) sign has to be considered for TO tangential(LO 
axial) modes. This expression, developed through a first- 
neighbor tight-binding model was then quantitatively con- 
firmed by direct DFT calculations 34 . The TO tangential/LO 
axial terms of the non-analytical part of the dynamical matrix 
can thus be written as: 



jjTO/LO _ 



■ E 



k dk ,n^(k>)] 



iu,i{k') - £v,f(k' + 
(11) 

where a factor 2 is included to take into account the contribu- 
tion of the two equivalent points K and K'. 

Since we are interested in phonons close to T, we should 
consider the q — > limit in Eq. [TT] 



B. Contribution of the v = bands 

The bands intersecting at the K point of the graphene BZ 
are at a good approximation linear around K, as shown in the 
previous section, with slopes ±/3 = 14.1 • ^ eV. For these 
branches, v = 0, so that the geometric contribution 1 ±008(6* + 
9') equals or 2. In practice, in the LO axial case it equals 
2 when the transition involves interband states with opposite 
slope, and vanishes for the intraband transitions; the opposite 
behavior occurs in the TO tangential case^. For LO axial 

/ + f^des Eq..[l2lbecomes , 

T i f $iqrT\£v. f ■ tu.i) cos(y + 9 )\ 



<1) 



£LO,0 



ngn\e u j ■ 



(G r >F 



dk' 



f{(3k') - f(-pk') 



(3k' 



In the limit T — * 0, one obtains 



D 



TO/LO 



4V3(G 



r/F 



dJfe' 



—k 




f[eu,i(k')}-f[e v ,f(k')] 



£LO,0 



4agV3 



i(k') ~ e vJ (k') 



[1 ± sign(e u j ■ e l 



i) cos(20)] 



(G r ) F ln 



\(3k\ 



(16) 



(17) 



E 



de 



vJ {k')] ■ [l±cos(20)] 



The quantitative result depends on the value of k, which en- 
ters the additive constant defined for lo lo at the end of the 
next subsection. Since the density 2>)f states is independent 
of energy for the v = band, the non-analytical direct EPC 
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contribution to the TO tangential mode is also a constant, in- 
dependent of e p 



£TO,0 



(18) 



which we include in the additive constant defined for loto- In 
other words, our model predicts, in the limit of low tempera- 
tures, a logarithmic divergence of the dynamical matrix at T 
for zero doping in the LO axial case, while the TO tangen- 
tial modes are simply red-shifted by the non-analytical EPC 
contribution. 



C. Contribution of the v = ±1 bands 

Let us now consider the second lowest v = ±1 bands. 
Since there is a large gap Eff between the ir and it* bands, 
the interband transitions cannot give rise to KAs, and are ex- 
pected to contribute negligibly to the Fermi level dependence 
of the phonon frequencies. On the other hand, in the case 
of doped nanotubes where the v = ±1 band is partly popu- 
lated, intraband transitions involve a density of states which 
diverges in one dimension, leading to possible Kohn anoma- 
lies. We develop the following analytical expressions based 
on the ideal hyperbolic v = ±1 bands. According to Eqs. 
|T2j[l3] and[14] the contribution of intraband transitions to -Dp 
amounts to 

^to/lci = _ 8 ^ pD0Si ^^ . [1 ± cos (26) F )] (19) 

where the v ± 1 degeneracy is taken into account by a factor 
of 2, and 



DOS 1> „.(e F )=0(e F -Ef 1 /2) 



n(3d 



lei - (E?{/2) 2 
(20) 

is the total DOS of the v = 1, ir* band (see Eq.O. The angular 
factor at the Fermi wavevector is 



cos(26>/ 



(21) 



where k[^, = ±-g <J ep — (Ef(/2) 2 is the Fermi momentum 
of the v = 1, 7r* band. We can notice that the EPC direct con- 
tribution D^P' 1 is proportional to the density of states of the 
v = ±1 conduction bands at the Fermi energy, modulated by 
the cosine factor. Kohn anomalies would thus occur where the 
density of states diverges. With a few algebraic calculations 
one obtains, for LO axial phonons, the compact form: 



D 



LO,l 



(GDf- 



(E?{/2y 



(e F -E%/2) 



(22) 



For TO tangential modes through similar calculations one fi- 
nally obtains: 



1J T 



o\/3 (G 2 



TT 2 d(i 



r/F- 



4-(#rY/2) 2 



(23) 



In this latter case, the divergence of the DOS is canceled by 
the cosine factor, which equals for ep = E\[j2. In the fol- 
lowing subsection, we will thus calculate through this model 
the r frequency of the TO/LO modes as: 



TO/LO 



nTO/LO,0 . nTO/LO,l 
/ TO/LO, other -.9 , + U T 

{uj r ) 2 + — 



M 



(24) 

where w TO/LO : othcr _ ^ ^to/lo con i; ams m e unper- 
turbed value uj-p = 1581 cm -1 , plus an additive constant 
C TO / LO , which also includes the terms discussed at the end 
of the previous subsection, and is chosen as to match the DFT 
phonon-EZF results at zero doping; its numerical value, for 
k = O.lfzt; i s 15.1/39.1 cm,- 1 in the TO/LO case respec- 
tively. 



D. Results: Kohn anomalies within the adiabatic 
approximations 

At this point, the actual determination of the effect of EPC 
in the vibrational properties of doped nanotubes can be in 
principle carried out by adopting approximated calculation 
schemes. The DFT phonon calculation of a SWNT is very 
demanding and is doable only at a relatively high electronic 
temperature, where the potentially interesting behaviors are 
smeared out. However, the comparison of the calculations 
done on a real SWNT with those done using the phonon- 
EZF (Fig. [T]) shows that a quantitative determination of the 
phonon frequency can be done neglecting the nanotube cur- 
vature. Now, we will analize an integral model which takes 
another step in identifying the key role played by graphene- 
derived EPC in the vibrational properties of SWNTs, and its 
low-temperature limit can be put into an analytical expression. 
We report in Fig.[3]the frequencies of the LO axial and the TO 
tangential modes of the G peak at an electronic temperature of 
315 K as a function of the Fermi energy, and thus of the elec- 
tronic doping, as calculated i) from the phonon-EZF method; 
ii) from the numerical integration of Eq.QT|at T=315 K (poly- 
nomial fits are used to describe the electronic bands); Hi) from 
its analytical limit at T=0 K, where we use the ideal conical 
bands (Eqs. El [H]|22l|23j and El, with E^f/2 red-shifted 
by 0.06 eV as explained in subsection lHEl 

A first inspection of the graphs indicates that the agreement 
within the different approaches, and in particular between the 
DFT calculations and the integral model, is excellent. The LO 
axial frequency drops at ep — due to a Kohn anomaly aris- 
ing from the first conduction band, and D^°'° varies as In cf 
according to Eq. [17] Another KA occurs when the Fermi en- 
ergy reaches the minimum of the second bands Ef[ /2 = 0.83 
eV. This second KA is much stronger than that at tp = 0, 
in agreement with Eq. [22] which predicts a variation go- 
ing as — 1 / Je 2 F — (E^/2) 2 . A frequency drop of about 40 

cm -1 , with respect to zero doping, was reported in a theoret- 
ical calculation of alkali-doped SWNTs 29 at a doping level of 
about N e i/C ~ 0.02, in good agreement with our results (see 
Fig- 03- Other KAs occur each time the Fermi level increases 
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FIG. 3: (Color online) Variation, within the adiabatic approximation, 
of the LO axial (black) and TO tangential (red) T frequencies with 
the Fermi level at T=315 K, calculated from DFT phonon-EZF (full 
diamonds), from the integral model (solid line, Eq. II It . and from 
the analytical model at T=0 K (dashed line, Eqs. [T7] [22] and 1231 . 
The metal/carbon composition corresponding to the Fermi level is 
indicated on the top of the figure. 



case, due to their intrinsic one-dimensional nature, and the 
use of the adiabatic BO approximation leads to significant de- 
viations from experimental results. This could be particularly 
true in the case of Raman scattering where a sinusoidal excita- 
tion induces the oscillatory motion of ions. In order to check 
whether a proper inclusion of non-adiabaticity would affect 
our previous results, we introduce a non-adiabatic model, 
based on time-dependent perturbation theory. In the process 
of absorption and emission of a phonon by the electrons, the 
phonon energy is no longer neglected: we replace, in Eq. QT| 
the energy difference between the two electronic scattering 
states,(e k ,i - £k+q,/)» bv ( £ k,i - £k+ q ,/ + ^ q + iS). Here 
we add a small imaginary part S to the energy to control the 
divergences. Eq. [TTIthen becomes: 



^TO/LO 



a 2 V3(G 



r/F 



E 



dk' 



-A- 



/K»(fc')]-/K/(fc' + g)] 
tv,i{k') - £v,f(k' + q) + huiq + it 

(25) 



The frequency is obtained as ujq = y - — j^-, that is, we take 
the principal part of the above integral by letting 5 — *■ 0. Since 
Dq depends on the frequency, this equation should be solved 
self-consistently. 

In the q — > limit, Eq.[T2]then becomes, in the NA case, 



such that a new conduction band is populated. In TO tangen- 
tial modes, due to the vanishing of the EPC interband matrix 
element, the first band does not give rise to a KA; the curve 
is practically flat until Ef[ /2. The effect of the second con- 
duction band induces a variation going as - J e 2 F — (Ef(/2) 2 
(Eq ~ 



However, as we will see in the following section, although 
the physics described by these results seems very intriguing, 
the main outcome of this comparison is the accuracy of the 
integral model, with respect to full DFT calculations, in pre- 
dicting the Raman properties of SWNTs, and in particular the 
anomalous doping dependence of the vibrational axial mode. 



IV. NON-ADIABATICITY 



^to/lo _ a 2 V3(Gl) F f k A1J J ^ / [e v ,i(k')] - f [e uJ {k')] 



D 



TT 2 d 



f c ]k' i V /t^ )]--* 



) + hu>r + i5 



(26) 

Here only the interband transitions contribute to the dynami- 
cal matrix, since, for an optical phonon, the limit 



lim 

q^o e v j(k') 



f[e»,f(k')]-f[e„Ak' + q)] 



u j(k' + q) + Tiujq + iS 



= 



(27) 



We can anticipate that the absence of the intraband contribu- 
tions dramatically affects the results when such terms are im- 
portant, i.e. when the Fermi level is close to a band 
minimum. As in the previous section, we study the contribu- 
tion of the v = 0, ±1 bands to the non-analytical part of the 
dynamical matrix. The mode frequencies are obtained from 
the dynamical matrix as in Eq.[24l 



A. Time-dependent perturbation theory 

The vast majority, if not the totality, of ab initio phonon- 
calculations rely on the so-called Born-Oppenheimer adia- 
batic approximation. In this framework one can decouple 
the electron motion from the ion dynamics, on the basis of 
their large mass and velocity difference, and thus treat the 
electronic properties as they were completely independent of 
the ionic motion. This approximation for phonon calculation 
is equivalent to first-order time-independent perturbation the- 
ory and is wholly justified for insulators and semiconductors. 
Even in ordinary metals a proper, specific treatment of the 
electronic degrees of freedom, such the inclusion of a finite 
Fermi-Dirac electronic temperature, usually suffices to avoid 
a failure of the BO approximation. 

Metallic SWNTs do however represent an "extraordinary" 



B. Contribution of the v = bands 

Eq.|26]can be integrated, by using linear bands and T = 0, 
analogously to the calculations developed in the static case; 
we thus obtain for the non-analytical part of Dr in the case of 
LO axial phonons 



TT 2 d(3 



In 



\2e F + hu\P\\2e F ~ hu\P\ 



\2f3k\* 



(28) 



while the TO tangential term TJ™' is zero since it involves 
the intraband terms only. We note that the KA observed in 
the static case for zero doping is replaced by two logarith- 
mic divergences at a doping level e F = ±fiuo^,° /2, which 
is N e i/C w 0.0015 (~ MC 65 o). As shown schematically in 
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FIG. 4: . Schematic representation of electronic transitions allowed 
by conservation of energy and momentum in electron-phonon scat- 
tering processes involving a phonon with momentum q and energy 
fcjq. An abrupt change in the phonon dispersion, i.e. a Kohn 
anomaly, occurs when, by change of the Fermi energy, one of this 
transitions becomes allowed or forbidden by the Pauli exclusion prin- 
ciple. In practice, this occurs when the Fermi level crosses a tip of 
one of the arrows corresponding to allowed transitions. Top panel: 
transitions at q=0 (r). Bottom panel: transitions at a finite q-point. 



the top panel of Fig. [4] T transitions conserving energy and 
momentum in the scattering process can indeed occur when 
connecting two bands separated by an energy equal to hur, 
which is in practice the case only in correspondence of the 
electronic vertical transitions (q=0) between the it and the it* 
bands. 



C. Contribution of the \ v I > 1 bands 



band terms are suppressed by the dynamical effects (Eq. 
in both LO and TO modes and thus /)LO,TO,t>^o _ q i n( j eec i 
in the hyperbolic band it is not possible to conserve both en- 
ergy and momentum in an electron-phonon scattering process 
with a q = optical phonon, see Fig. |4] In addition to the 
analytical calculation at T=0, we performed numerical inte- 
grations of Eq.|26]at different temperatures (77 K, 315 K). As 
in the static case, we used polynomial fits to describe the band 
dispersions. We report in Fig. [5]the variation of the LO axial 
mode as function of the Fermi level in the non-adiabatic case 
at those two temperatures. The results are consistent with the 
analytical ones obtained at vanishing temperature; the only 
Kohn anomaly is observed, as explained in the previous sub- 
section, at ?kJr/2. All the other Kohn anomalies, that in the 
adiabatic approximation and thus in all standard DFT calcu- 
lations are due to the population of new energy levels, com- 
pletely disappear. This is a general behavior for vertical tran- 
sitions. 



D. LO phonon linewidth 

An important measurable quantity is the phonon linewidth. 
In ref.— it is shown that it can be split into a EPC di- 
rect contribution, only relevant in metallic nanotubes, and a 
term due to inhomogeneous broadening and anharmonicites, 
common to all nanotubes and estimated around 10 cm -1 
from experimental data. The LO axial phonon linewidth 
(full width at half maximum) j^q ' as function of the elec- 
tronic doping can instead be calculated 39 from the imagi- 
nary part of the non-adiabatic dynamical matrix, Eq. [25] as 
7q = \^(Dq)\/ (wqM). An identical result is obtained by 
using the Fermi golden rule, as in refill. The result of ref. 24 
derived in nanotubes for zero doping can be easily generalized 
for any doping as: 



7LO i e F) 



TrdfiMuj^ 
79 cm nr, 



(G 



T/F 



1 



e (e F -hw^°/2)/k B T + 1 
1 



e (ep- hu%°/2)/k B T + 1 e (e F +h^°/2)/k E 



At zero doping this FWHM equals, at T=315 K, about 60 
cm- 1 for a (9,9) SWNT. We show the EPC LO linewidth in 
the top panel of Fig. [5] At T — ► this term is constant up 
to the Kohn anomaly at a N el /C « 0.0015 (MC 650 ) doping 
level, and vanishes abruptly for higher doping levels, while 
a smoother behavior is observed at finite temperature. Since 
all the terms contributing to the EPC part of the dynamical 
matrix vanish for the TO tangential modes, the corresponding 
phonon linewidth will vanish as well. 



The effect of the non-adiabatic terms in the v = bands 
contribution merely consists in a splitting of the zero-doping 
Kohn anomaly observed in the static case. The picture 
changes completely for the \u\ > 1 bands. In fact, given the 
energy gap between the tt and tt* bands having v ^ 0, only the 
intraband terms are relevant to Kohn anomalies. These intra- 



E. Finite-q results 

In the previous subsections, we have developed our analyt- 
ical model at the T point of the Brillouin zone, and we have 
shown that the proper treatment of non-adiabatic terms lifts 
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FIG. 5: (Color online) Lower panel: variation of the LO axial (G~ 
peak in Raman experiments, black) and TO tangential (G + peak, red) 
r frequencies with the Fermi level at T=315 K (solid line) and at 77 
K (LO axial, dashed line), calculated from the non-adiabatic integral 
model (Eq.l25t. Top panel: variation, at the same temperatures, of 
the EPC contribution to the FWHM of the LO axial mode with the 
doping level. The metal/carbon composition corresponding to the 
Fermi level is indicated on the top of the figure. 



all the divergences resulting from intraband, vertical transi- 
tions. However, as shown schematically in the lower panel of 
Fig. |U at finite q the outcome is different from the T case. 
Indeed, non-vertical intraband transitions conserving energy 
and momentum can occur when q is approximately larger than 
oj q /vM, vm being the maximum electronic band velocity. We 
report in Fig. [6] the dependence on electronic doping of the 
LO axial frequency at q=0.05 • — and 0.1 • — , calculated nu- 

^ J jl . an an 

merically through Eq. |25]at T=77 K and 315 K. In this case 
the intraband Kohn anomalies observed in the adiabatic ap- 
proximation are not suppressed, as it is the case for vertical 
transitions. The interband KA is shifted towards higher dop- 
ing levels with respect to the T case. The intraband KAs are 
shifted in the same direction with respect to the energy levels 
Both interband and intraband Kohn anomalies appear 
in doublets, separated by hujq. 

By using the expressions for EPC derived in Refsi2i2& for 
the scattering involving phonons near the K point, it can be 
shown that the behavior of the highest optical phonon at a 
wavevector q+K is very similar to the one reported in Fig. [6] 
for the LO axial frequency at q. This phonon near the K point 
is responsible for the D band observed experimentally in Ra- 
man spectra. 
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FIG. 6: (Color online) Lower panel: variation, with respect to zero 
doping, of the LO axial frequencies with the Fermi level at T=315 K 
(black) and T=77 K (red), calculated from the non-adiabatic integral 



model (Eq.|25t, The phonon momentum is equal to q=0. 1 



Top 



panel: same, at q=0.05 ~. The dotted vertical line indicates the 
minimum of the second conduction band. The energy values /xi-, 
fJ,2-, and jU2+ correspond to the levels schematically shown in 
Fig-El A similar behavior is expected for the highest optical phonon 
near the K point, which is responsible for the Raman D band. The 
metal/carbon composition corresponding to the Fermi level is indi- 
cated on the top of the figure. 



V. DISCUSSION 
A. Physical properties of metallic nanotubes 

In this work we presented a combination of ab initio DFT 
calculations with an integral and analytical model developed 
to take into account in a simple way the direct EPC contribu- 
tion to the dynamical matrix, responsible for the divergences 
originating the KAs. Since however the mathematics behind 
our models might be cumbersome, here we summarize the 
most interesting physical properties of metallic SWNTs we 
observe, as well as the quantitative results on experimentally 
measurable quantities. 

We showed that KAs in SWNTs can be identified by con- 
sidering electron-phonon scattering processes that conserve 
both energy and momentum (see Fig. [4] and relative discus- 
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sion). Since there are strong disagreements between the adi- 
abatic and the non-adiabatic case, only the results obtained 
through the latter model are reliable. The failure of the adia- 
batic calculation originates from the fact that, in this approx- 
imation, the phonon energy is neglected and the energy con- 
servation in the electron-phonon scattering processes is thus 
violated. Since the curves shown in Figs.Q]and[3]are obtained 
through the adiabatic model and through DFT frozen-phonon 
and DFPT calculations, those results should not be used for 
comparison to experimental data. On the other hand, all the 
results presented in Figs.[5]and [6] which conflict with any ex- 
isting DFT calculation based on the BO approximation, are in 
principle correct, and are thus meant to be directly compared 
to experiments. 

The results shown in Fig.[5]promise to be useful tools to de- 
termine experimentally the doping level of individual metal- 
lic SWNTs. In particular, by measuring the shift of the G~ 
peak with respect to zero doping and its linewidth, one can 
determine the variation of the Fermi energy and thus the dop- 
ing level. Those results refer to a (9,9) nanotube of diameter 
1.23 nra; however, they can be easily generalized to any d- 
diameter metallic SWNT by observing that the dependence of 
both the linewidth of the G~ peak and its frequency shift with 
respect to zero doping case are proportional to 1 jd and inde- 
pendent of the chirality. We note here that we are assuming 
that the G peak originates from single resonance. 

It was showed that carbon nanotubes 36 ' 37 , as graphit o 40 ' 41 , 
present a free-electron-like interlayer state whose energy level 
decreases rapidly with increasing doping, and reaches the 
Fermi energy at about the saturation doping of ^MCs- We 
showed that KAs occur at T only when the Fermi energy 
is in proximity of a band crossing with non-zero interband 
electron-phonon coupling. Besides the crossing of the ir- 
bands at zero doping, we estimate that the first band-crossing 
with symmetry-allowed transitions is that between the inter- 
layer level and the ir* (v = 0) bands. Such crossing is a 
candidate for a Kohn anomaly. This could explain the strong 
softening of about 50 cmT 1 observed at saturation doping in 
experimental Raman studies on the G peak in alkali-doped 
SWNTs^. 

Finally, we predicted KAs to occur in phonons with finite 
momentum q not only in proximity of a band crossing, but 
also each time a new band is populated (see Fig. [6]). Therefore, 
by experimentally probing such KAs, the filling of the hyper- 
bolic bands could be detected. Phonons at finite-q vectors are 
experimentally accessible to double-resonant Raman scatter- 
ingi 3 -, and correspond in the Raman spectra, e.g., to the defect- 
activated D peak at about 1350 cm -1 , and to the second- 
order (two-phonon) 2D and 2G peaks at approximately 2700 
and 3200 cmT 1 . We showed in Fig. [6] the variation of the 
phonon frequency with respect to doping at q=0.05 • |^ and 
0.1 ■ These phonon momenta are experimentally accessi- 
ble to double-resonant Raman scattering 19 . 

Although our computational study has been limited to a 



metallic SWNT, our analytical results allow to predict the be- 
havior of the high-energy Raman peaks in semiconducting 
nanotubes. Since in this case only intraband terms contribute 
to the EPC term of the dynamical matrix, we expect the G 
peak to be unaffected by electron doping. On the other hand, 
analogously to the metallic case, the defect-activated D peak 
should be sensitive to doping, provided that the excitation en- 
ergy is sufficient to populate the first excited electronic level 
in the conduction band. 



B. Conclusions 

In this work we presented a theoretical study on the vi- 
brational/Raman properties of electron-doped SWNTs us- 
ing Density-Functional Theory and analytical models. We 
performed our calculations within the adiabatic Born- 
Oppenheimer approximation, but we also used time- 
dependent perturbation theory to explore non-adiabatic ef- 
fects beyond this approximation. We showed that the Born- 
Oppenheimer approximation, predicts, for increasing doping 
levels, a series of EPC-induced KAs in the vibrational mode 
parallel to the tube axis at the T point of the Brillouin zone, 
usually indicated in Raman spectroscopy as the G~ peak. 
Such Kohn anomalies would arise each time a new conduc- 
tion band is populated. However, we showed that they are 
an artifact of the adiabatic approximation, which is the stan- 
dard approach for ab-initio phonon calculations. The inclu- 
sion of non-adiabatic effects dramatically affects the results, 
predicting KAs at T only when ep is close to a band cross- 
ing Ex- For each band crossing a double KA occurs for 
€f = Ex ± hui/2, where hui is the phonon energy. In par- 
ticular, for a 1.2 nm metallic nanotube, we predicted a KA 
to occur in the so-called G~ peak at a doping level of about 
N el /C = ±0.0015 atom (e F « ±0.1 eV) and, possibly, 
close to the saturation doping level (N e i/C^ 0.125), where 
an interlayer band crosses the ir* nanotube bands. Further- 
more, we predicted that the Raman linewidth of the G~ peak 
significantly decreases for \ep | > foui/2. Thus our results 
provide a tool to determine experimentally the doping level 
from the value of the frequency shift and from the linewidth 
of the G~ peak. Finally, we predicted Kohn anomalies to 
occur in phonons with finite momentum q not only in prox- 
imity of a band crossing, but also each time a new band is 
populated. Such Kohn anomalies should be observable in the 
double-resonant Raman peaks, such as the defect-activated D 
and peak, and the second-order peaks 2D and 2G. We also 
predict that in semiconducting nanotubes the G peak should 
be insensitive to doping, while the D peak should be affected 
analogously to the metallic case. 

We thank N. Bendiab, S. Piscanec and A. C. Ferrari for 
useful discussions. DFT calculations have been performed at 
the IDRIS French National Computational Facility under the 
projetCP9-61387. 
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